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Abstract 


A design strategy for optimal design of com- 
posite grid-stiffened panels subjected to global 
and local buckling constraints is developed using 
a discrete optimizer. An improved smeared stiff- 
ener theory is used for the global buckling anal- 
ysis. Local buckling of skin segments is assessed 
using a Rayleigh-Ritz method that accounts for 
material anisotropy and transverse shear flexibil- 
ity. The local buckling of stiffener segments is 
also assessed. Design variables are the axial and 
transverse stiffener spacing, stiffener height and 
thickness, skin laminate, and stiffening configura- 
tion. The design optimization process is adapted 
to identify the lightest-weight stiffening config- 
uration and pattern for grid-stiffened composite 
panels given the overall panel dimensions, design 
in-plane loads, material properties, and bound- 
ary conditions of the grid-stiffened panel. 
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Nomenclature 

a Axial stiffener spacing 

b Transverse stiffener spacing 

h Stiffener height 

t Skin laminate thickness 

t s Stiffener thickness 

LAM I Design variable for stacking 

sequence of skin laminate 
ICON Design variable for stiffening 
configuration 

Ac? Global buckling load factor 

X s k Buckling load factor for skin 

segment 

A x Buckling load factor for axial 

stiffener segment 
A 2 Buckling load factor for 

transverse stiffener segment 
A3 Buckling load factor for diagonal 

stiffener segment 

Introduction 

An aircraft in flight is subjected to air loads 
associated with all flight conditions including 
maneuver and gust conditions. These external 
loads are resisted by the structure, and an in- 
ternal load distribution is established based on 
the structural layout and external loads. These 
internal loads, which depend on their location 
in the aircraft structure, may cause either over- 
all panel buckling of stiffened panels, buckling 
of the skin between stiffeners, or stiffener crip- 
pling. Hence, an efficient and accurate method 
for developing a buckling-resistant design of gen- 
eral stiffened composite panels is needed to iden- 
tify the most effective grid-stiffened geometries 
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for structural panels subjected to combined in- 
plane loading conditions at different locations in 
fuselage and wing structures. The identification 
of structurally efficient grid-stiffened geometries 
also requires integration of optimization tech- 
niques with accurate structural analyses. Op- 
timization of composite structures has been of 
considerable interest in recent years. Composite 
plates have been optimized to maximize buckling 
loads and a summary of this work is presented in 
the literature survey of Reference 1 . According to 
Reference 1 , the basic design problem is to deter- 
mine the stacking sequence of a composite lam- 
inate since, in many practical applications, the 
ply orientations are limited to either 0, 90 or ±45 
degrees and, the laminate thicknesses can only 
be integer multiples of commercially available ply 
thicknesses. Thus, the optimization of a laminate 
stacking sequence involves discrete design vari- 
ables associated with manufacturing constraints 
and represents an integer programming problem. 

The optimum design of stiffened panels that 
satisfies buckling constraints has also been of in- 
terest in recent years (e.g., Refs. 1-6). These 
studies did not consider stiffener spacing as de- 
sign variables even though skin thickness, stiff- 
ener thickness and stiffener height were used as 
design variables. For the most part, gradient- 
based optimizers were used in References 3-6. 
However in Reference 2, the ranking method was 
used as the optimizer, while in Reference 1, a ge- 
netic algorithm 7,8 was used to optimize the lam- 
inate stacking sequence in the skin and stiffener 
elements. Geodesically stiffened panels were con- 
sidered in References 3 and 6, while orthogrid 
panels were considered in Reference 2. Axially 
stiffened panels were considered in References 1 
and 5. The optimization of grid-stiffened panels 
with stiffener spacing and stiffener layout as dis- 
crete design variables poses the same problems 
as that of optimizing the laminate stacking se- 
quence discussed in Reference 1 and significantly 
expands the design space. Stiffening configura- 
tion and stiffener spacings are discrete variables 
since the optimum grid-stiffened geometry may 
contain any combination of axial, transverse and 
diagonal stiffeners, the stiffener spacing can only 
occur as certain multiples of the length and width 
of the panel, and the stiffener and skin thick- 
nesses can only be integer multiples of commer- 
cially available ply thicknesses. The genetic al- 
gorithm has emerged as a viable tool for deal- 
ing with the problem of discrete variables and 
with the need to find multiple minima. A ge- 


netic algorithm evolves the design by randomly 
searching the design space and maintaining a 
family of designs for each generation (or itera- 
tion) of the design. This process provides multi- 
ple near-optimum designs for evaluation and se- 
lection rather than a single-design configuration 
provided by gradient-based algorithms. 

The present paper presents the analysis and 
design strategy for grid-stiffened composite pan- 
els subjected to combined loads and a global 
buckling design constraint. The global buckling 
constraint for the grid-stiffened panel is impor- 
tant to prevent localized skin buckling. This con- 
straint reduces the loss of aerodynamic perfor- 
mance caused by the buckling of wing or fuse- 
lage surfaces and prevent the failure of a panel 
by stiffener-skin separation 9 after buckling. 

Panel Buckling Analysis 

The analysis and design of grid-stiffened 
composite panels subjected to combined loads 
require several key steps. In the present study, 
acceptable designs are those which buckle glob- 
ally and do not exhibit any local skin buckling 
or stiffener crippling. The first step in the design 
process is to assess the global buckling response 
of a grid-stiffened panel. Once this global buck- 
ling response is determined, the second step is 
to determine the local skin buckling response for 
the quadrilateral and/or triangular skin segments 
between the stiffeners. The third step is to deter- 
mine whether stiffener buckling or stiffener crip- 
pling has occurred at this global buckling load 
level. This sequence of steps is performed re- 
peatedly in a design cycle until an optimum or 
near-optimum design is obtained. 

The global buckling analysis is based on a 
Rayleigh-Ritz method using a first-order, shear- 
deformation theory and an improved smeared- 
stiffener modeling approach discussed in Refer- 
ence 10. The buckling analysis of local skin- 
segments is also based on a Rayleigh-Ritz analy- 
sis using a first-order, shear-deformation theory 
and accounts for material anisotropy. Bound- 
ary restraints on the skin segments are pro- 
vided by the stiffeners and hence, the analy- 
sis must be capable of accommodating a va- 
riety of boundary conditions and a variety of 
skin segment shapes. 11,12 In most cases, the skin 
segments for grid- stiffened panels will have ei- 
ther a general parallelogram-shaped or a general 
triangular-shaped planform. 

In addition to analyzing the local skin seg- 
ment for buckling, the local stiffener segments 
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must be analyzed to determine whether stiffener 
crippling will occur, 4 Accordingly, the stiffener 
segment at the nodes or intersection points of 
the stiffeners are assumed to be clamped while 
the stiffener-skin attachment is assumed to be a 
simple support. 

These global and local analysis methods 
have been integrated into a computer code to pro- 
vide a computationally efficient tool for predict- 
ing the buckling load of grid-stiffened composite 
panels. 

Panel Design Procedure 

The design of grid-stiffened composite pan- 
els requires that many of the design variables, 
such as stiffener spacing and stiffener thicknesses, 
may only have certain discrete values rather than 
varying continuously over design space. Also, a 
“family” of acceptable designs is often needed for 
selecting designs that accommodate manufactur- 
ing constraints rather than a single-point design. 
Gradient-based methods for structural optimiza- 
tion are not appropriate in this case, since they 
lead to a single-point design. 

Over the last several years, researchers have 
investigated the use of genetic algorithms as a 
method for “evolving” a given design problem to 
a family of near-optimum designs (e.g., see Ref- 
erences 1, 7 and 8). Based on Darwin’s theory of 
the survival-of-the-fittest , the genetic algorithm 
involves the random creation of a design popu- 
lation that “evolves” towards some definition of 
fitness. The genetic algorithm is attractive due to 
the simplicity of its approach using discrete vari- 
able combinatorics. The genetic algorithm can be 
used directly to solve unconstrained optimization 
problems, while constrained optimization prob- 
lems must first be transformed into an uncon- 
strained optimization problem (e.g., use of an ex- 
terior penalty function). Stochastic processes are 
used to generate an initial population of individ- 
ual designs, and the algorithm then applies the 
principles of natural selection and survival of the 
fittest to find improved designs. Furthermore, 
since the discrete design procedure works with a 
population of designs, it can explore a large area 
of design space and identify multiple minima or 
maxima. This attribute is a major advantage 
since the converged solution contains many op- 
tima of comparable performance. The cost of 
having a large number of function evaluations 
is offset by the fact that many near-optimum 
designs are now available. In a gradient-based 


optimization procedure, only a single-point de- 
sign, usually the extremum closest to the start- 
ing point, is obtained. However, different start- 
ing points can be tried to increase the chance of 
locating the global optimum as well as other lo- 
cal optima. The genetic algorithm produces a 
population or family of good designs which may 
include the global optimal design, rather than a 
single design. The above mentioned advantage 
of the genetic algorithm makes it an appropriate 
optimizer that can be exploited in developing a 
design optimization tool for general grid-stiffened 
panels. 

Design Problem Definition 

The present design problem is to minimize 
the weight per unit area of a grid-stiffened com- 
posite panel given the design loading condition, 
the length and width of the panel, the material 
properties for the skin and stiffeners, and the 
boundary conditions of the panel. As shown in 
Figure 1, a general grid-stiffened panel may be 
considered as an assembly of a repetitive unit or 
a unit cell. A unit cell is a repetitive unit of 
skin and stiffener elements. The design variables 
include stiffener spacings (a, 6), skin stacking se- 
quence, stiffening configuration, stiffener thick- 
ness (< 4 ), and stiffener height (h\ = h 2 = h 3 = 
h). The stiffener spacings (a, 6), stiffener height 
(/i), and stiffener thickness (t s ) are shown in the 
figure of the unit cell of Figure 1. Also shown in 
Figure 1 is the skin thickness (/) which depends 
on the skin stacking sequence. The stiffening con- 
figuration depends on the combination of axial, 
transverse and diagonal stiffeners in the unit cell 
(Figure 1). For manufacturing and assembly rea- 
sons, all stiffeners are assumed to be of the same 
height and thickness. The design sought here is 
a panel of minimum weight in a certain design 
space which buckles globally at the design loads 
This design problem can be defined by setting up 
the optimization procedures in the following way. 
First, the global buckling load is assumed to be 
a scalar multiple of the design loads and has the 
form 

(N X )cr = (Ny) cr = A(^A 7 y, 

( N X y) cr = ^G^xy ( 1 ) 

where N Xi N yi N xy are the applied in-plane pre- 
buckling loads. These values represent the design 
loads for the grid-stiffened panel. Second, the 
design constraints imposed on the panel include 
requirements that 
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1. The critical buckling load should be 
greater than or equal to the design loads, 
that is, Xg > 1. 

2. The skin segments should not buckle at 
the critical buckling load, that is, X sk > 1. 

3. The stiffener segments should not crip- 
ple at the critical buckling load, that is, 

A 2 , A 3 > 1 where Ai, A 2 , A 3 are 
the crippling load factors of the x-direction 
stiffener, y-direction stiffener and diagonal 
stiffener, respectively. 

The general form of each constraint equation is 
written as 

9j — “ 1) < 0.0 j = l,...,N c (2) 

Finally, the “Fitness” expression based on an ex- 
terior penalty function approach is 


Fitness = ( 


Max- 


Q 


Q 

F(x } n) 


)= 


^(X) + r,Ef [|ffj(X)| + #(X)P 


( 3 ) 


where 

X = design variable vector 
F(X y n) = modified objective function 
W(X) = weight of panel per unit area 

r « [|5;( x )l + 5j( x )] 2 = penalty 

function 

Q = normalizing constant 

N c = number of design constraints 

n = penalty parameter 

i = generation or iteration cycle in the 

optimization procedure. 

Once the global buckling load factor has 
been determined, the loads acting on the stiff- 
ener and skin segments have to be determined by 
distributing the loads among the skin and stiff- 
eners based on their extensional stiffnesses. The 
procedure for distributing the applied loads for a 
general grid-stiffened panel is discussed in Refer- 
ences 6 and 13. 

The weight per unit area of the grid-stiffened 
panel is 


w = ^(tui + w 2 + W 3 + w 3 ) 

where 


wi — 2 h a t s 

w 2 — 2 hbt s 

W3 — 2 h t s \fa 2 + 6 2 

w s = a b t (4) 

wi is the volume of the axial stiffeners in the unit 
cell, w 2 is the volume of the transverse stiffeners 
in the unit cell, iu 3 is the volume of the diagonal 
stiffeners in the unit cell, w s is the volume of the 
skin in the unit cell, t is the thickness of skin, and 
p is the mass density of the material. 

Design Process Based 
on a Genetic Algorithm 

Implementation of the genetic algorithm is 
shown schematically in Figure 2. The design pro- 
cess begins with a random selection of a specified 
number of designs which comprise the initial pop- 
ulation (i.e., first generation) for the genetic al- 
gorithm. The problem parameters such as mate- 
rial properties, length and width, boundary con- 
ditions of the panel, and design loads are input 
into the analysis processor routine. The buckling 
analyses are performed which provides the crit- 
ical eigenvalues for the global buckling response 
of the grid-stiffened panel, and the local buckling 
response of the skin and stiffener segments. The 
weight per unit area of the grid-stiffened panel 
is also computed. This procedure is repeated for 
each design configuration in the population. The 
“fitness” processor then evaluates the “fitness” of 
each design using Equation (3) and assigns a rank 
based on the fitness expression or the objective 
function. The current population of design con- 
figurations is then processed by the genetic oper- 
ators (crossover, mutation, and permutation) to 
create a new population of design configurations 
for the subsequent generations which combines 
the most desirable characteristics of the previ- 
ous generations. Designs from previous genera- 
tions may be replaced by new ones (i.e., children) 
except for the “most fit” designs (i.e., parents) 
which are always included in the next generation. 
The process is repeated until design convergence 
is obtained, which is defined herein by specifying 
a maximum number of generations that may oc- 
cur without any improvement in the best design. 
The design procedure is demonstrated on flat and 
curved grid-stiffened panels in the following sec- 
tions. 

Numerical Results for Flat 
Grid-stiffened Panels 

A 20.0-in.-long and 56.0-in.-wide flat grid- 
stiffened composite panel representative of a 
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generic transport helicopter fuselage structural 
component is designed to demonstrate the capa- 
bilities of the present design optimization tool. 
The panel is optimized for a load case with 400.0 
lbs/in. of axial compression, which is considered 
to be the critical load case for this panel. The 
termination criterion is 30 generations, and the 
population size is eight. The probabilities used 
for crossover, mutation, and permutation are 1.0, 
0.10, and 0.95, respectively. Also, the penalty 
parameter r* is kept constant for all iterations 
since the genetic algorithms maximize Equation 
(3) more efficiently. The skin laminate stacking 
sequence selected is [±45/90/0] a . The stiffeners 
are made of 0° material only. The nominal ply 
mechanical properties used are: E\\ = 20.2 Msi; 
E 22 = 1.9 Msi; G 12 = G \3 = G 23 — 0.73 Msi and 
V 12 = 0.3. The mass density of the material p, is 
0.0570 lbs/in. 3 The grid-stiffened, panel has sim- 
ply supported boundary conditions on all edges, 
and the skin segments are also considered to be 
simply supported. The grid-stiffened panel is as- 
sumed to have only axial and diagonal stiffeners. 
The axial and transverse stiffener spacings con- 
sidered are such that all stiffener patterns closely 
approximate an isogrid configuration. Hence, the 
axial stiffener spacing a and transverse stiffener 
spacing b are not independent but are consid- 
ered as a single design variable, (i.e., (a 3 ,6 s ) is 
one design variable). The stiffener height h and 
thickness t s are also design variables. The design 
space explored is indicated in Table 1, where the 
height of the triangle (b s / 2) is kept between 2.9 
and 6.0 in., and the stiffener aspect ratio ( h/t s ) 
is kept between 4.5 and 9.0 due to manufacturing 
constraints. Each design variable is permitted to 
assume only eight discrete values by the genetic 
algorithm FORTRAN code being used. 

The five best designs from the optimization 
are shown in Table 2. These designs buckle glob- 
ally at the indicated load factor of A <3, since the 
local buckling load factors (A 3 j &, Ai, A3) are all 
greater than one. The first and third designs 
have a global buckling load factor Xq of 0.995 
and 0.991, respectively, and could still represent 
acceptable designs. The second and third designs 
have different stiffener spacing, and yet the global 
buckling load factor and weight per unit area of 
these two panels are very close to one another. 
The fourth and fifth designs also, have similar re- 
sults. Hence, the ability of the genetic algorithm 
to obtain multiple optima of comparable perfor- 
mance is demonstrated. The best design for the 
grid-stiffened panel with a [±45/90/0] 3 skin lam- 


inate is the first design given in Table 2. Design 
convergence was achieved after 40 iterations. 

The grid-stiffened panel is now assessed us- 
ing a representative flight-load conditions given 
in Table 3. The results are also shown in Ta- 
ble 3. Of the six load cases shown in Table 3, 
load cases 2 to 5 are the flight-load cases. The 
first and sixth load cases are used to complete the 
shear and axial-load interaction curve. The panel 
buckles globally for most of the load cases con- 
sidered. Crippling of the diagonal stiffener occurs 
for the load case of N xy = 100 Ibs/in. and N x 
= 10 lbs/in. In addition, the grid-stiffened panel 
exhibits local skin buckling for the load case of 
N x = 174.0 lbs/in. and N xy = 154.0 lbs/in. The 
global buckling load factor Xq for this load case 
is 2.2382, and the skin local buckling load factor 
A,* is 0.9855. The buckling load factor is 2.2057 
for the load case of N x = 174.0 lbs/in. and N xy 
= 154.0 lbs/in. 

Numerical Results for Curved 
Grid-stiffened Panels 

The fuselage design of a generic wide-body 
transport aircraft is typically divided into four 
different quadrants. These quadrants include a 
crown panel, two side panels and a keel panel. 
A side-quadrant panel is considered herein and 
designed for global buckling. The side-quadrant 
panel is shown in Figure 3 and longerons and 
frames divide the side-quadrant panel into four 
curved panels. Each panel is chosen to be 22.0- 
in. long and 22.0-in. wide with a radius of 120.0 
in. in the width direction. The dimension of 22.0 
in corresponds to the frame spacing. Panel 1 is 
the forward-top panel of the side-quadrant panel 
and is subjected to N x = 1250 lbs/in., N xy = 250 
lbs/in., and N y = -2200 lbs/in. (hoop tension). 
Panel 2 is the aft-top panel of the side-quadrant 
panel and is subjected to N x = 300 lbs/in., N xy 
= 1350 lbs/in., and N y = -2200 ibs/in. (hoop 
tension). Panel 3 is the bottom-top panel of the 
side-quadrant panel and is subjected to N x = 
2250 lbs/in., N xy = 250 Ibs/in., and N y = - 
2200 lbs/in. (hoop tension). The panel hoop 
tension is due to internal pressurization of the 
fuselage. The nominal ply mechanical proper- 
ties used are: E\\ = 20.2 Msi; E 22 = 1-9 Msi; 
G \2 — G13 = G23 = 0-73 Msi and v \2 = 0.3. The 
mass density of the, material p is 0.0570 lbs/in. 3 
The grid-stiffened panel is assumed to have sim- 
ply supported boundary conditions on all edges, 
and the individual skin segments are considered 
to be simply supported also. The stiffeners are 
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made of unidirectional material. The termination 
criterion for the design evolution is 25 generations 
with no improvement in the “best ” design, and 
the population size is twelve. The probabilities 
used for crossover, mutation, and permutation 
are 1.0, 0.10, and 0.95, respectively. Also, the 
penalty parameter r, is kept constant for all iter- 
ations. 

The design variables are the axial stiffener 
spacing (a), the transverse stiffener spacing (6), 
the stiffener height (ft), the stiffener thickness 
(t s ), the stacking sequence of the skin lami- 
nate (LAM I), and the stiffening configuration 
(ICON) which is a design variable indicating 
the combination of axial, transverse, and diag- 
onal stiffeners in a unit cell. The design space 
explored for a, 6, ft, and t s is shown in Table 
4 for Panels 1 and 3, and in Table 5 for Panel 
2. The design space for LAM I and ICON is 
described in Table 6 for Panels 1, 2 and 3, e.g., 
when LAM I = 1, then the skin stacking sequence 
is [±45/0] 2 5 and when ICON =1, then the stiff- 
ening configuration consists of axial stiffeners. In 
either design space, the minimum stiffener spac- 
ing is restricted to two inches, and the aspect 
ratio of the stiffener (h/t s ) is kept between 3.5 
and 10.5 due to manufacturing constraints. Each 
design variable can assume eight discrete values. 
A modified Sanders-Koiter shell theory is used to 
account for transverse shear deformation in the 
Rayleigh-Ritz buckling analysis. 13 

The results obtained for Panel 2 using the 
present optimization tool are shown in Table 7. 
The panel designs presented in Table 7 buckle 
globally at the corresponding global load factor 
of \q. The genetic algorithm produces a large 
pool of acceptable designs in this case. Most 
of the acceptable designs only have axial stiffen- 
ers. These axially stiffened panels have stiffener 
spacings similar to those of the first three de- 
signs presented in Table 7, but with different val- 
ues of stiffener height and thickness. The fourth 
design in Table 7 has transverse and diagonal 
stiffeners (ICON — 6); however, this design is 
14 percent heavier than the first design. Panels 
stiffened in multiple directions have redundant 
load paths and typically exhibit better damage 
tolerance characteristics than panel stiffened in 
one directions. Therefore, the design process was 
repeated using a modified design space for the 
stiffening configuration where values of ICON 
= 1 and 2 are replaced by ICON — 5 and 6, 
respectively. That is, all designs include stiffen- 


ers in multiple directions. The results for this 
optimization are shown in Table 8. The pan- 
els presented in Table 8 buckle globally at the 
corresponding global load factor of A G . The ge- 
netic algorithm produces a large pool of accept- 
able designs with axial and transverse stiffeners 
(ICON = 3). These panels have stiffener spac- 
ings represented by the first four designs of Table 
8, with variation in stiffener height and thickness 
and with a skin laminate stacking sequence of 
[±45/0] 2 , (LAM I = 1) or [±45/90] 2 , (LAMI 
= 2). The weight of each panel is comparable 
to the weight of the axially stiffened panels pre- 
sented in Table 7. The fifth design has both axial 
and diagonal stiffeners (ICON = 5), while the 
sixth design is the same as the fourth design pre- 
sented in Table 7. The best design for Panel s is 
the design with transverse and diagonal stiffen- 
ers, which is the fourth design in Table 7 and the 
sixth design in Table 8. This design is preferred 
over the lighter weight axially stiffened panel and 
the axially and transversely stiffened panel since 
the panel with diagonal stiffening may be more 
damage tolerant. 9 The sixth design in Table 8 is 
also preferred over the fifth design in Table 8, 
since it has more stiffeners than the axially and 
diagonally stiffened panel and its weight is close 
to that of the latter. The convergence behavior 
of the genetic algorithm for Panel 2 is shown in 
Figure 4. The convergence is faster when there 
is no restriction on the stiffening configuration 
variable (ICON). The optimization process is 
performed for Panels 1 and 3, and the best de- 
sign obtained for these two panels are shown in 
Table 9. 

The results for the best design obtained for 
Panel 3 for the design load cases used for Panel 1 
and Panel 2 are shown in Table 10. For the loads 
corresponding to Panel 1, the panel buckles glob- 
ally at a buckling load factor of 2.114, and hence, 
this design represents a conservative design for 
these loads. For the loads corresponding to Panel 
2, the buckling load factor of the diagonal stiff- 
ener is 0.255, and, hence, the buckling deforma- 
tion contains local buckling of the diagonal stiff- 
ener at a load factor of A 3 x X G — 1.172. There- 
fore, this design also represents a conservative de- 
sign for the loads for Panel 3. Hence, the best 
design for Panel 3 can also be used for Panel 1 
and Panel 2 with a weight penalty of 6.7 and 22.9 
percent, respectively, when compared to their re- 
spective best design. For a grid-stiffened fuse- 
lage structure construction, it is desirable that 
all the panels have the same stiffener pattern due 


6 



to both manufacturing and structural considera- 
tions. The region where stiffeners from two dif- 
ferent panels connect (heavy lines in Figure 3) 
have to be joined and is a critical area especially 
if the stiffening patterns of the two panels are 
not the same. To alleviate that problem, these 
regions have to be reinforced, which could result 
in a higher weight penalty. 

Concluding Remarks 

A minimum-weight design optimization tool 
with buckling constraints has been developed 
for grid-stiffened panels using global and local 
buckling analyses and a genetic algorithm. De- 
sign variables used are axial and transverse stiff- 
ener spacing, stiffener height and thickness, skin- 
laminate stacking sequence, and stiffening config- 
uration. Results for flat and curved grid-stiffened 
panels indicate that the best design configura- 
tions obtained by the genetic algorithm depend 
mostly on the load cases. For most combined ap- 
plied load cases, the best designs have diagonal 
stiffeners. However in certain cases, the best de- 
sign obtained by the genetic algorithm, based on 
buckling considerations only, may not be suitable 
for structural applications due to manufacturing, 
joining, and damage-tolerance considerations. In 
this case, the pool of acceptable designs obtained 
by the genetic algorithm provides alternative de- 
sign options for a given structural application. 
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Table 1 Design space for flat grid-stiffened pan- 
els. 


Stiffener 
spacing, in. 
( a a > b s ) 

Stiffener 
height, in. 

(h) 

Stiffener 
thickness, in. 
(*.) 

(6.667 , 11.200) 

0.49375 

0.060 

(5.714 , 10.182) 

0.50000 

0.066 

(5.000 , 8.615) 

0.50625 

0.072 

(4.444 , 8.000) 

0.51250 

0.078 

(4.444 , 7.467) 

0.51875 

0.084 

(4.000 , 7.000) 

0.52500 

0.090 

(3.636 , 6.222) 

0.53125 

0.096 

(3.333 , 5.894) 

0.53750 

0.102 


Table 2 Best designs obtained by the genetic 
algorithm for a grid-stiffened panel. 


Design 

variables 

Buckling load 
factors 

weight, 

lbs/ft 2 

a 9 = 5.000 in., 
b s — 8.615 in., 
h = 0.5375 in., 
t s = 0.06 in. 

A a = 0.995, 
A sfc = 1.024, 
A! = 1.341, 
A 2 = 21.208 

0.578 

a s = 4.444 in., 
b s = 8.000 in., 
h = 0.5375 in., 
t s = 0.06 in. 

A g = 1.047, 
A s * = 1.179, 
A, = 1.470, 
A 2 = 25.608 

0.596 

a s = 4.444 in., 
b s = 7.467 in., 
h = 0.5125 in., 
t s — 0.06 in. 

A g = 0.991, 
A sk = 1.392, 
Ai = 1.637, 
A 2 = 24.507 

0.594 

a s = 4.444 in., 
b s = 7.467 in., 
h = 0.5125 in., 
t s = 0.072 in. 

A g = 1.139, 
A , k = 1-285, 
Ai = 2.146, 
A 2 = 32.120 

0.634 

a $ = 4.444 in., 
b s = 8.000 in., 
h = 0.5375 in., 
i s = 0.072 in. 

A g = 1.207, 
A, fc = 1.084, 
Ai = 1.920, 
A 2 = 33.458 

0.636 


Table 3 Buckling loads for grid-stiffened panels. 

Load case 
number 

N x , 

lbs/in. 

Nxy j 
lbs/in. 

Buckling load 
factor, A 

i 

400.0 

0.000 

0.995“ 

2 

396.0 

130.0 

1.002“ 

3 

326.0 

178.0 

1.211“ 

4 

271.0 

139.0 

1.458“ 

5 

174.0 

154.0 

2.206 i> 

6 

10.00 

100.0 

4.145“ 


a global buckling. 

6 skin buckling. 
c diagonal stiffener buckling. 
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Table 4 Design space for a, b, h and t, of Curved 
Panels 1 and 3. 


a, in. 

6, in. 

/i, in. 

t,, in. 

11.000 

11.000 

0.49375 

0.060 

8.800 

8.800 

0.50000 

0.066 

7.333 

7.333 

0.50625 

0.072 

6.286 

6.286 

0.51250 

0.078 

5.500 

5.500 

0.51875 

0.084 

4.889 

4.889 

0.52500 

0.090 

4.400 

4.400 

0.53125 

0.096 

4.000 

4.000 

0.53750 

0.102 


Table 5 Design space for a, b, h and t a of Curved 
Panel 2. 


a, in. 

6, in. 

h f in. 

t 3 , in. 

11.000 

11.000 

0.30 

0.042 

8.800 

8.800 

0.32 

0.048 

7.333 

7.333 

0.34 

0.054 

6.286 

6.286 

0.36 

0.060 

5.500 

5.500 

0.38 

0.066 

4.889 

4.889 

0.40 

0.072 

4.400 

4.400 

0.42 

0.078 

4.000 

4.000 

0.44 

0.084 


Table 7 Best designs obtained by the genetic 
algorithm for grid-stiffened Curved Panel 2. 


Design 

variables 

Weight, 

lbs/ft 2 

Buckling load 
factors 

b = 11.0 in., 
h = 0.360 in., 
t s = 0.060 in., 
LAMI = 
[±45/0] 2 „ 
ICON = 1 

0.6232 

A <3 = 1.061, 
X,k = 16.494 
Aj = 2.745 

b = 7.333 in., 
h = 0.340 in., 
t , = 0.042 in., 
LAMI = 
[±45/0] 2 „ 
ICON = 1 

0.6229 

A g = 0.998 
A jfc = 221.17 

Aj = 1.62 

6 = 6.2857 in., 
h = 0.340 in., 
t s = 0.042 in., 
LAMI = 
[±45/0] 2 „ 
ICON = 1 

0.6283 

A g = 1.129 
A sk = 206.79 
Aj = 1.45 

a = 5.5 in., 
b = 6.2857 in, 
h = 0.320 in, 
t s = 0.054 in, 
LAMI = 
[±45/0] 2 „ 
ICON = 6 

0.711 

A g = 0.996, 
A , k = 5.01, 
A 2 = -.30, 

A 3 = 1.650 


Table 6 Description of design space for ICON 
and LAM I for Panels 1, 2, and 3. 


Integer 

value 

LAMI 

ICON 

i 

[±45/0] 2j 

axial stiffeners 

2 

[±45/90] 2j 

transverse stiffeners 

3 

[±45/0/90] 2j 

axial and transverse 
stiffeners 

4 

[±45/0 2 ] 2j 

diagonal stiffeners 

5 

[±45/90 2 ] 2 , 

axial and diagonal 
stiffeners 

6 

[±45/0 2 /90] 2j 

transverse and 
diagonal stiffeners 

7 

[±45/0/90 2 ] 2s 

axial, transverse and 
diagonal stiffeners 

8 

[±45/0 2 /90 2 ] 2j 

no stiffeners 
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Table 8 Best designs obtained by the genetic 
algorithm for grid-stiffened Curved Panel 2, with 
ICON = 1, 2 excluded from the design space. 


Design 

variables 

Weight 

lbs/ft 2 

Buckling load 
factors 

a = 1L0 in., 

6 = 8.8 in., 
h = 0.340 in., 
t s = 0.042 in., 
LAM I = 
[±45/0] 2j| 
ICON = 3 

0.6389 

A g = 1.038, 
Ajjt — 3.05, 

Aj = 2.01, 

A 2 = -0.149 

a = 11.0 in., 

6 = 11.0 in., 
h = 0.340 in., 
t s = 0.048 in., 
LAM I = 
[±45/90] 2s , 
ICON = 3 

0.6396 

A g = 1.022, 
Ajjt = 2.07 
A t = 1.28, 
A 2 = -0.360 

a = 11.0 in., 

6 = 7.333 in., 
h = 0.320 in., 
t s = 0.042 in., 
LAM I = 
[±45/90] 2 , , 
ICON = 3 

0.6411 

A G = 1.076, 
Ajjt = 4.141 
Xi = 1.07 
A 2 = -0.366 

a = 8.8 in., 

6 = 8.8 in., 
h = 0.340 in., 
t s = 0.042 in., 
LAM I = 
[±45/90] 2 ,, 
ICON = 3 

0.6442 

A g = 1.083, 
X sk = 2.979 

Aj = 1.08, 

A 2 = -0.300 

a = 11.0 in., 
b = 11.0 in., 
h = 0.340 in., 
t s = 0.084 in., 
LAM I = 
[±45/0] 2 „ 
ICON = 5 

0.6938 

A g = 2.160, 
X sk = 2.442 
Ax = 3.23, 
A 3 = 1.234 

a = 5.5 in., 
b = 6.2857 in, 
h = 0.320 in, 
t s = 0.054 in, 
LAM I = 
[±45/0] 2j , 
ICON = 6 

0.711 

A g = 0.999, 
X ak = 5.01 
A 2 = -.30 
A 3 = 1.650 


Table 9 Best designs obtained by the genetic 
algorithm for grid-stiffened Curved Panels 1 and 
3. 


Design Weight, Buckling load 

variables lbs/ft 2 factors 

Panel 1 

a = 4.8890 in., 0.819 A G = 1.080, 

b = 6.2857 in., X sk = 1.290, 

h = 0.500 in., A 2 = -0.394, 

t, = 0.060 in., A 3 = 18.189 

LAM I = 

[±45/90] 2j 

ICON = 6 

Panel 3 

a = 5.500 in., 0.874 A G = 1.014, 

b = 4.400 in., X, k = 1.369, 

h = 0.50625 in., A 2 = -0.781, 

t, = 0.072 in., A 3 = 1.191 

LAM I = 

[±45/90] 2j , 

ICON = 6 


Table 10 Buckling loads for the best design of 
Panel 3 subjected to the load cases of Panels 1 
and 2. 

Panel Load, Buckling load 

number lbs/in. factors 


Panel 1° N x = 1250, A G = 2.114, 

N v = -2200, X, k = 1.488, 

N xy = 250. A 2 = -0.374, 

A 3 = 1.0113, 

Panel 2 b N x = 300, A G = 4.595, 

N y = -2200, X sk = 2.409, 

N xy = 1350. A 2 = -0.172, 

A 3 = 0.255 

“Global buckling at A G = 2.114 
‘Local buckling of diagonal stiffener 
St A critical — 1.172 
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Figure 1: Unit cell of a grid-stiffened panel show- 
ing design variables. 




j 


Figure 2: Flow chart for the optimization using 
the genetic algorithm. 


Figure 3: Side quadrant panel of fuselage struc- 
ture. 




Figure 4: Design convergence for composite grid- 
stiffened curved panel (Panel 2). 
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Figure 1 Unit cell of grid-stiffened panel 
showing design variables. 




Figure 3 Side quadrant panel of a fuselage 
structure. 


Figure 2 Row chart for the optimization 
procedure using the genetic algorithm. 
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Figure 4 Design convergence history for a 
composite grid-stiffened curved panel (Panel 2). 


